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1.  Introduction 


A  petroleum  product  called  fog  oil  is  the  primary  substance  examined.  Fog 
oil  is  best  characterized  as  a  light  oil  simUar  to  that  used  to  lubricate  small 
machinery.  For  many  years,  fog  oil  has  been  atomized  to  produce  a  fog¬ 
like  aerosol,  most  commonly  seen  in  skywriting  but  also  used  for  many 
years  as  an  Army  obscurant,  smoke. 

The  purpose  of  this  study  was  to  determine  absorption  and  scattering  of 
electromagnetic  radiation  by  this  type  of  hydrocarbon  aerosol  in  an 
atmospheric  transmission  window  at  wavelengths  between  3  and  4  pm. 
Although  interest  is  shown  in  describing  interaction  of  fog  oil  aerosol  with 
high-energy  laser  beams,  the  optical  coefficients  are  available  only  in 
unpublished  reports  and  the  values  at  the  wavelengths  of  interest  immensely 
vary.  In  this  report,  the  values  of  n  and  k  for  the  midinihued  (IR)  region 
are  shown  to  be  in  reasonable  agreement  with  those  of  the  most  recent 
unpublished  report. 

The  measurements  of  optical  properties  are  made  on  fog  oil  aerosols 
dispersed  and  continuously  mixed  in  controlled  environmental  chambers. 
Aerosol  mass  densities  and  size  distributions  are  determined  by  using 
dosimetry  and  commercially  available  techniques.  Settling  characteristics 
of  the  well-mixed  chamber  aerosols  yield  supportive  information  regarding 
the  particle  size  distributions.  Mie-theory  ^culations  for  the  measured 
size  distributions  are  based  on  values  of  the  complex  refractive  indices 
measured  by  Weng.  [1]  The  study  extends  to  optical  properties  for  other 
similar  petroleum  products. 
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2.  Approach 


The  ai^roach  consisted  of  two  stages; 

1.  A  combination  photoacoustical  and  extinction  measurement  was 
performed  at  the  single  wavelength  of  3.39  /xm.  Fog  oil  was  nebulized 
in  a  controlled  environment  (chamber)  where  the  absorption  and 
extinction  of  the  aerosol  could  be  determined  in  situ.  Aerosol  density 
and  size  distribution  were  concurrently  measured.  From  the  optical  and 
dosimetry  measurements,  the  mass  normalized  absorption  and  total 
scattering  coefficients  were  calculated. 

2.  A  relative  extinction  profile  was  measured  from  2.5  to  approximately 
12  ^m  by  an  IR  scanning  transmissometer  (IRST).  The  aerosol  density 
and  size  distribution  were  determined. 

The  Mie  calculations  of  the  extinction  efficiencies  at  3.39  /xm  were  based 
on  the  previously  measured  aerosol  size  distributions  and  the  referenced 
complex  indices.  Agreement  with  the  photoacoustic  portion  of  the 
experiment  was  good  and  justified  the  method  of  normalizing  the  relative 
extinction  profile.  The  result  is  a  continuous,  absolute,  extinction  profile 
of  the  aerosol  over  the  region  of  interest. 

A  series  of  Mie  calculations  spanning  the  larger  range  of  2.5  to  12  /xm, 
using  the  measured  size  distribution  and  referenced  indices  of  fog  oil, 
provide  comparison  with  the  IRST  results. 
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3.  Environmental  Chamber  for  Measurements  at  3.4  /an 


Fog  oil  was  atomized  in  a  small  volume  (0.07  m^)  chamber  by  a  De  Vilbis 
pharmaceutical  nebulizer.  A  small  fan  established  and  maintained  a  nearly 
uniform  distribution  of  particles  throughout  the  chamber.  The  stirring 
ensured  maintenance  of  the  spatial  distribution  but  was  not  vigorous  enough 
to  significantly  alter  the  size  distribution.  The  photoacoustic  system 
continuously  sampled  ft'om  this  plenum. 
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4.  Dosimetry  and  Size  Characterization 


Dosimetric  measurements  were  made  throughout  each  experiment  to  record 
the  aerosol  density  as  a  function  of  time.  A  Gelman  filter  holder  and  filter 
were  inserted  into  a  port  projecting  into  the  chamber.  Gelman  type 
AE  fiberglass  filters  with  a  O.l-fim  pore  size  were  used.  The  volume  flow 
rate  (approximately  3  L/min)  and  the  sample  time  were  recorded.  The  flow 
rate  and  sample  times  were  chosen  to  minimize  disturbances  and  allow 
collection  of  a  measurable  mass  amount.  The  filters  were  weighed  on  a 
Metder  precision  balance.  Typical  mass  samples  were  several  milligrams 
with  repeatability  of  0.2  mg. 

The  size  distribution  for  the  chamber  was  characterized  in  a  separate  set  of 
measurements.  Two  types  of  particle  spectrometers  were  used  to 
characterize  the  distribution:  (1)  Particle  Measuring  System  (PMS)  model 
ASASP-X  (0.07  <  radius  <  3  /xm,  with  small  orifice)  and  (2)  CSASP-1(X) 
(0.2  <  radius  <  16  fitn).  A  typical  size  distribution  measured  in  the  fixed 
frequency  portion  of  this  experiment  is  shown  in  figure  1.  The  number 
density  distribution  with  radius  is  very  strongly  peaked  near  a  radius  of 
0.1  fxm,  and  few  particles  exceed  a  radius  of  5.0  urn. 
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5.  Photoacoustic  Measurements 


The  photoacoustic  and  extinction  measurements  were  made  with  a  flow¬ 
through  photoacoustic  unit,  described  previously  by  Bruce  and  Pinnick  [2] 
and  Bruce  and  Richardson,  [3]  containing  a  3.39  /xm  HeNe  laser  source, 
beam  positioning  optics  and  a  modulator,  a  cylindrical  acoustical  cavity 
with  an  isolating  array  of  acoustic  filters,  and  a  detector. 

The  laser  beam  is  directed  down  the  axial  center  of  the  cavity  by  two 
positioning  mirrors  (figure  2).  The  laser  beam  is  modulated  by  a 
mechanical  tuning  fork  chopper  at  the  fundamental  longitudinal  resonant 
frequency  of  the  acoustic  cavity.  The  signal  from  the  microphone 
embedded  centrally  in  the  acoustic  cavity  wall  is  amplified  with  phase- 
sensitive  circuitry.  Kreutzer  [4]  and  Trusty  [5]  show  that  this  signal  is 
linearlv  prqportional  to  the  amount  of  energy  absorbed  by  the  aerosol.  The 
consta..i  proportionality  is  repeatedly  determined  by  a  calibration  gas 
(isopropyl  alcohol  for  3.4  /tm),  for  which  the  only  significant  attenuation 
is  due  to  absorption. 


WINDOW 


Figure  2.  Gas  and  aerosol  spectrophone/transmissometer  for  0.63, 1.1,  or  3.4  fim. 

The  beam  intensity  is  monitored  by  a  detector  at  the  opposite  end  of  the 
cavity.  The  attenuation  by  the  aerosol  is  used  to  calculate  the  extinction 
coefficient  according  to  Beer’s  law,  with  a  forward  scattering  correction 
quantified  by  Deepak  and  Box.  [6] 
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Aerosol  flow  through  the  photoacoustical  cavity  is  regulated  and  varies 
from  1  to  3  m/s.  During  measurement,  a  valve  system  alternately  selected 
fresh  air  and  aerosol  flow  through  the  spectrophone.  This  was 
systematically  performed  to  monitor  the  unattenuat^  laser  power  and 
establish  an  intensity  baseline.  After  a  uniform  flow  was  achieved,  the 
laser  power  and  spectrophone  signal  were  recorded  as  a  function  of  time. 

From  the  photoacoustic  portion  of  the  experiment  the  absorption  and 
extinction  coefficients  as  a  function  of  time  (figure  3)  were  measured.  The 
shaded  regions  represent  the  dosimetric  sampling  time  intervals.  The  log 
averages  were  taken  within  the  intervals  to  determine  an  associated 
absorption  and  extinction  coefficient  (via  the  absorption  calibration 
constant).  An  aerosol  density  was  computed  for  each  interval  on  the  basis 
of  the  mass  samples  obtained.  Mass  normalized  efficiencies  were 
calculated  from  the  absorption  and  extinction  coefficients  and  the  aerosol 
densities  at  specific  times.  Hie  forward  scattering  correction  for  the 
measurement  geometry  and  aerosol  size  distribution  of  figure  1  is  not 
significant. 


Figure  3.  Time  evolution  of  extinction  and  absorption  coefficients. 


Tlie  average  measured  and  calculated  efficiencies  at  3.39  /xm  are  shown  in 
table  1.  The  calculated  values  are  obtained  using  equation  (1)  and  the 
measured  density  distribution. 
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where 

£  =  extinction  coefficient 

p  =  aerosol  density 

—  number  of  particles  per  radius  increment 

i 

Qi  =  Mie  coefficient 

Pb  —  bulk  density  of  the  aerosol  particle 


Table  1.  Fog  oil  attenuation  efficienci^  at  3.39  iim  for  the  size  distribution  in 
figure  1 


Measured 

Calculated 

Absoiption 

0.483  ±0.03 

0.401  ±0.06  (mVg) 

Total  Scattering 

0.887  ±0.05 

0.809  ±0.12  (mVg) 

The  uncertainties  given  in  table  1  are  the  one-sigma  variations.  The 
uncertainties  quoted  for  the  calculated  values  are  solely  derived  from 
variations  in  the  size  distribution.  The  measured  absorption  coefficient  is 
17  percent  higher  than  calculated,  and  the  total  scattering  coefficient  is 
9  percent  higher  than  calculated.  Both  are  within  the  respective  propagated 
uncertainties,  but  it  is  likely  (and  is  illustrated  later  in  this  report)  that  a 
significant  portion  of  the  error  results  from  an  uncertainty  in  the  imaginary 
component  of  the  index  of  refraction  at  3.39  pm. 
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6.  Spectrally  Continuous  Measurements 


Spectrally  continuous  measurements  were  conducted  in  a  larger  (cubic 
meter)  chamber  with  two  large  fans  (operating  at  controlled,  reduced  rates 
of  rotation)  producing  an  aerosol  circulation  pattern.  Two  axial  tubes, 
opposite  each  other  across  the  center  of  the  chamber,  acted  as  the 
transmission  windows.  Fog  oil  was  nebulized  in  the  same  manner  as  the 
photoacoustic  measurements  in  an  effort  to  produce  a  size  distribution  close 
to  that  of  the  0.07-m^  chamber.  Three  jdiarmaceutical  nebulizers  were  used 
to  atomize  the  fog  oil  and  produced  a  particle  size  distribution  similar  to  the 
photoacoustic  measurement  (figure  4).  Because  the  size  distributions  for 
both  chambers  indicate  that  nearly  all  the  mass  was  attributable  to  particles 
with  radii  above  0.1  fim  (as  illustrated  in  figure  7),  the  measured  mass  did 
not  require  correction  for  filter  pore  size.. 

A  1(X)0-K  glow  bar  was  used  as  an  IR  source.  Radiation  was  collimated 
and  focused  through  the  chamber  and  onto  the  IRST. 

The  IRST  utilizes  a  continuously  variable  circular  filter  whose  resolution 
depends  on  the  wavelength  and  the  filter  quadrant.  The  wavelength 
resolution,  generally  adequate  for  aerosol  spectra,  was  just  barely  adequate 
for  this  study  at  1.3  percent  of  the  3.39-/tim  wavelength. 

The  filtered  IR  radiation  was  detected  by  a  Mercury-Cadmium-Telluride 
detector.  The  signal  was  amplified  (lOOOX)  and  processed  by  an 
8085  microprocessor,  and  noise  reduction  schemes  such  as  boxcar 
integration  were  applied. 


17 


0-0  o.j  1.0  10.0 


RADIUS  (//ml 

Figure  4.  F(^  oil  particle  size  distributions  for  the  two  chambers. 


A  relative  extinction  curve  was  computed  from  2.S  to  4.0  fan  using  data 
measured  by  the  IRST.  The  curve  was  normalized  to  the  measured 
absolute  extinction  efficiency  at  3.39  ^m. 

Figure  S  shows  a  comparison  of  the  calculated  and  measured  (mass 
normalized)  extinction  coefficients  for  the  large  chamber. 

18 


l» 

u 

z 

ly 


Figures.  Comparison  of  calculated  (Mie  tlieory)  and  experimeotal  results. 


TTic  line  laser  absorption  measurement  fells  within  the  strong  adsorption 
band  indicated  in  figure  5.  Structurally,  this  band  results  primarily  from 
symmetric  and  asymmetric  stretch  bands  of  CHj  and  CHj.  Spectra 
presented  in  Szymanski  [7]  show  that  the  spectra  of  fog  oil  are 
characteristic  of  many  hydrocaibons  and  other  substances.  The  net  foim 
is  presumed  to  be  a  superposition  of  the  contributions  caused  by  the  liquids 
of  the  constituent  molecules  weighted  by  their  relative  occurrence  in  the 
presence  of  the  net  liquid.  Quantitatively,  the  problem  is  even  more 
complex  than  a  superposition  of  the  constituents. 

It  is  presumed  that  discrepancies  between  the  calculated  and  measur^ 
curves  arose  from  uncertainties  in  the  measured  size  distribution,  noise  in 
the  spectral  measurements,  uncertainties  in  the  index  measurements,  and 
chemical  differences  between  materials  affecting  the  indi^.  Fog  oil  is  a 
term  that  describes  a  range  of  qjecifications  for  liquid  hydrocaibons 
identified  in  a  report  by  Katz  et  al.;  [8]  therefore,  actual  indices  may  vary. 

The  previously  discussed  Mie  calculations  were  expanded  to  include  the 
absorption  efficiencies  and  were  extended  out  to  12  pm  in  the  plots  in 
figure  6.  The  plots  show  that,  for  the  aerosol  size  distribution  produced 
(probably  typical  for  fog  oil),  absorption  is  a  minor  contribution  to  the 
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ABSOLUTE  MIE  CALCULATION  EFFICIENCIES 


extinction  in  the  IR,  The  absofpdon  lesonance  in  the  vicinity  of  3.4  /zm  for 
which  the  absorption  is  ai^noximately  one-third  of  dw  total  extinction 
provides  the  largest  contribution. 


Figure  6.  Calculated  extinction  and  absorption  spectra  of  fog  oil. 
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7.  Effective  Radius  of  Fog  Oil  Particles 


A  rather  useful  and  independent  calculation  was  made  to  corroborate  the 
experimental  results.  The  calculation  involves  a  method  of  determining  the 
radii  of  particles  most  representative  relative  to  the  extinction  cross  section. 
The  approach  is  based  on  an  analysis  described  by  Bruce  et  al.  [9]  of  the 
settling  properties  by  the  time  dependent  decay  of  the  extinction  coefficient. 
The  method  considers  the  rate  at  which  particles  of  a  given  size  settle  out 
of  a  well-mixed  medium,  which  can  be  related  to  the  aerosol  aerodynamic 
cross  section.  The  result  yields  the  radius  with  the  largest  contribution  to 
the  associated  cross  section  (and  dirough  the  form  of  time  dependence  can 
give  a  measure  of  the  size  distribution  breadth).  Settling  should  not  be 
confused  with  diffusion,  a  much  slower  process  for  this  size  distribution. 
The  effective  radius  was  determined  using  equation  (2). 


C) 


where 

c  =  ®  (*  Jt  V) 

S  =  slope  of  the  optical  coefficient  versus  time 

h  =  settling  height  of  the  chamber 

V  =  viscosity  of  air 

g  =  acceleration  due  to  gravity. 

This  calculation  of  the  most  representative  radius  was  made  for  the 
photoacoustic  and  the  IRST  measurements  (because  the  size  distributions 
were  slightly  different)  and  the  f  values  for  the  extinction  cross  sections 
were  found  to  be  1.58  ±0.2  fim  ^hotoacoustic,  0.07-m^  chamber)  and 
2. 10  ±0.2  fim  (IRST,  1 .0-m’  chamber). 

The  f  values  were  compared  with  the  calculated  cross  sections  for  each 
case.  The  calculated  radii  correspond  very  well  to  the  respective  curve 
maxima  and  further  support  the  integration  of  the  two  experiments 
(figure  7). 
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Figure  7.  Relative  extinction  cross  section  as  a  function  of  radius  for  fog  oil 
aerosol. 
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8.  Inversion  of  the  Mie  Process  to  Determine  Complex 
Index  at  3.39  laa 


Although  the  intent  was  to  compare  measured  valiMS  of  absorption  and  total 
scattering  with  calculated  vidiKs  from  the  size  distribution  and  die 
previously  measured  complex  index,  it  would  perh^  lake  more  sense  to 
reverse  the  process  and  accurately  determine  the  values  of  n  and  k  by  an 
iterative  i^roach  using  the  Mie  theory.  The  inverse  reproach  has  been 
used  for  several  aeroMls  at  various  wavelengths  and  found  it  to  stabilize  fen* 
the  real  and  imaginary  components  when  a  substance  exhibits  significant 
absorption  and  scattering. 


Figure  8  illustrates  the  character  of  the  inversion  process  for  fog  oil  at 
3.39  im  using  f  as  the  effective  particle  radius.  Selection  of  a  unique  pair 
of  values  is  apparent  for  the  range  of  indices.  The  components  determined 
are  compared  with  the  previously  measured  values  in  table  2.  Use  of  the 
complete  size  distribution  provides  a  more  accurate  determination,  but  the 
example  in  table  2  illustrates  the  uniqueness  of  the  result. 


0.00  0.05  0.10  .123 


Scailering 


Absorption 


0.15 


Imaginary  Index 

figure  8.  Determination  of  complex  index  using  Nfie  theory,  effective  radius, 
measured  absorption,  and  total  scattering  coeflients  for  fog  oil.  Nfie  theory 
coefficients  form  an  increasing  progression  as  real  index  ircreases  (for  the  family  of 
corves)  in  increments  of  0.15. 
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TaUe  2.  Complex  index  using  the  absorption  and  scattering  ^fidencies,  dfective 
radius,  and  Mie  theory  as  compared  wi^  those  of  Weng.  [1] 


N 

K 

Calculated 

1.425 

0.125 

Weng  [1] 

1.450 

0.091 
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9.  Chemical  and  Optical  Analysis  of  Similar  Substances 


Fog  oil  and  several  similar  petroleum  products  were  analyzed  and  an 
attempt  was  made  to  relate  fog  oil  optical  properties  to  tlw  optical 
properties  of  other  liquids  possessing  similar  hydrocarbon  structures. 

A  Fourier  transform  IR  spectrometer  determined  the  bulk  absorption 
profiles  of  liquid  fog  oil,  common  diesel  fuel,  and  kerosene.  The  three 
substances  exhibited  nearly  identical  absorption  structure  within  the  IR 
(figure  9a). 

The  chemical  analysis  of  generic  petroleum  products  is  complex.  Even 
relatively  simple  oils  may  be  composed  of  many  species.  In  the  present 
investigation,  it  was  sufficient  for  comparison  to  characterize  each 
substance  by  its  overall  mass  profile.  Gas  chromatography  with  mass 
spectroscopy  was  employed  for  this  procedure.  Both  the  diesel  fuel  and 
kerosene  had  average  mass  peaks  between  the  C,2  and  hydrocarbons, 
and  the  less  volatile  fog  oil  processed  a  heavier  average  mass  peak  that 
ranged  from  to  C24  (figure  9b). 

An  excellent,  thorough  chemical  characterization  was  conducted  on  several 
different  varieties  of  fog  oil  by  Katz  et  al.  [8]  A  brief  summary  of  their 
results  follows: 

The  results  show  that  fog  oil  consists  of  nearly  pure  hydrocarbons,  with 
the  predominant  structures  being  mixtures  of  aliphatic  and  aromatic 
components  in  almost  equal  amounts.  Also  detected  were  small 
amounts  of  alcohols,  organic  acids,  and  esters  with  very  small  traces  of 
organic  nitrogen  derivatives,  llie  aliphatic  hydrocarbons  were  in  the 
C12  to  C22  range  and  the  aromatics  consisted  of  1-  through  4-member 
rings,  also  within  the  same  range.  Ibe  results  also  show  that  all  fog 
oils  contained  traces  of  copper  and  zinc,  with  copper  near  40  ppb  and 
zinc  varying  between  20  and  100  ppb.  Densities  ranged  between  0.89 
and  0.93  g/ml. 

A  distinctive  color  change  from  clear,  light  yellow  to  dark  brown  over 
a  24-month  period  was  noted. 
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Transmissivity  (a)  and  mass  spectra  (b)  of  liquid  fog  oil,  kerosene,  and  diesel  fuel 


10.  Comments  on  the  Size  Distribution 


Measurements  were  carefully  taken  to  to  ensure  that,  while  the  medium  was 
well  mixed,  the  size  distribution  was  minimally  perturbed  in  the  sampling 
process. 

The  instrumental  response  function  of  the  PMS  light  scattering  counters  to 
the  radii  of  the  spherical  particles  is  multivalued  in  the  resonance  region; 
therfore,  several  different  sized  particles  can  yield  the  same  response. 
According  to  Pinnick  and  Auvermann,  [10]  the  multivalued  response  can 
create  an  artificial  ripple  structure  in  the  size  distribution  unless  particles 
are  categorized  in  radius  ranges  designed  to  avoid  causing  an  artificial 
ripple  structure.  The  reduction  in  resolution  resulting  from  the  adjustment 
of  radius  size  increments  was  not  severe  for  the  measured  distributions  of 
fog  oU. 
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11.  Summary 


The  absorption  and  total  scattering  efficiencies  for  fog  oil  aerosol  were 
determined  and  compared  with  calculated  values  for  given  size 
distributions.  After  various  parametric  studies  and  analyses,  it  is 
considered  that  most  of  the  uncertainty  resides  in  the  size  distributions  and 
the  indices  (the  calculated  values). 

The  optical  measurements  were  extended  to  a  continuous  wavelength  span 
in  the  vicinity  of  3.39  nm  by  measuring  extinction  spectra.  Absorption 
spectra  were  not  measured  because  available  continuous  sources  do  not 
have  sufficient  power  to  drive  the  photoacoustical  measurement. 
Normalization  of  the  spectrum  using  the  laser  result  gains  credibility 
because  the  calculated  vdues  are  in  reasonable  agreement  with  the  form  of 
the  extinction  spectrum  and  with  the  absorption  and  extinction  values  at  the 
laser  line  wavelength. 

Indicating  similarity  of  optical  properties  for  related  liquid  hydrocarbons 
was  of  interest  and  was  investigated  here  only  for  the  bulk  materials,  but 
the  absorption  spectra  near  3.39  fim  were  found  to  be  similar  for  the  three 
substances  (liquid  fog  oil,  common  diesel  fuel,  and  kerosene). 

The  measure  d  size  distribution  was,  to  a  degree,  substantiated  by  agreement 
with  the  independent  settling  theory  calculations.  The  calculation  of  the 
effective  radius  for  an  optical  coefficient  agreed  well  with  the  peak  radius 
for  each  distribution  and  indicated  the  narrow  breadth  of  the  distribution 
through  the  small  curvature  of  the  semilog  plot  showing  decay  of  the  optical 
coefficient  with  time. 

The  effective  radius  and  the  coefficients  provided  a  basis  for  an  independent 
and  accurate  determination  of  the  complex  index  at  the  line  laser 
wavelength. 
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